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ABSTRACT  (Continue  on  ravaraa  alda  If  nacaaaary  and  Identify  by  block  numbor) 


This  report  describes  the  capabilities,  limitations,  and  operation  of  the  OSU, 
GTD-AMP  Computer  Code,  which  has  been  developed  to  calculate  the  electromagnetic 
radiation  from  currents  on  thin  wires  in  the  presence  of  a perfectly-conducting 
polygonal  plate  positioned  over  a reflecting  surface  of  Infinite  extent. 

Although  the  code  is  based  on  the  Geometrical  Theory  of  Diffraction  (GTD/  it 
is  designed  to  be  compatible  with  the  Antenna  Modeling  Program  (AMP).  The  OTD- 
AMP  program  described  herein  is  capable  of  handling  many  antenna-scattering 
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problems  once  the  currents  on  the  wire radiators  are  known. Ordinarily  these 
currents  are  found  from  the  AMP  program  in  the  manner  described  in  the  text; 
however,  they  may  be  specified  or  obtained  by  other  means.  This  combination  of 
the  moment  method  and  GTD  makes  it  possible  to  calculate  radiation  patterns 
which  could  not  be  calculated  by  either  method  individually.. 
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EVALUATION 


This  contractual  effort  resulted  in  the  development  of  a geometrical 
theory  of  diffraction  (GTD)  computer  code  using  the  output  of  the  antenna 
modeling  program  (AMP).  This  new  GTD-AMP  program  analytically  predicts 
the  far  field  pattern  of  wire  antennas  above  finite  size  ground  planes 
taking  into  account  the  mutual  coupling  between  current  elements  on 
the  wire  structure  and  the  ground  plane.  The  far  field  pattern  also  takes 
into  account  the  edge  diffraction  from  the  finite  ground  plane. 

The  results  of  this  effort  provide  a means  of  simulating  the  antenna 
and  its  environment  more  accurately  than  previously  accomplished  in  the 
past.  This  fits  into  the  RADC  Technology  Plan  (TPO-5)  for  analytically 
simulating  antenna  systems  to  determine  their  effectiveness  before  fabrication 
and  test. 

The  GTD-AMP  computer  code  is  presently  operational  on  the  RADC-HIS 
635  computer  facility  and  will  be  used  to  analyze  the  performance  characteris- 
tics of  radar  and  communication  type  antenna  systems. 
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I.  PROGRAM  CAPABILITIES 

This  document  describes  a computer  program  (referred  to  hereafter  as 
the  OSU  or  GTD-AMP  program)  which  has  been  developed  for  the  computation  of 
the  far  field  radiation  from  an  arbitrary  system  of  short  current  carrying 
wire  segments  situated  above  a flat  polygonal  plate  which  itself  is  located 
above  a planar  impedance  surface  or  a dielectric  medium  bounded  by  a plane 
surface.  The  wire  radiator-plate  geometry  is  shown  in  Figure  1. 

The  wire  radiator  may  have  up  to  300  segments  (this  number  could  easily 
be . increased) , which  can  be  arbitrarily  connected.  The  complex  currents 
(e^wt  time  convention)  flowing  on  each  wire  segment  are  assumed  to  be  already 
known.  The  wire  segments  may  be  arbitrarily  located  over  the  plate  provided 
they  are  at  least  X/4  away  from  the  nearest  edge. 

We  have  shown  that  if  the  wire  segments  are  located  over  the  surface  of 
the  plate  so  that  it  is  possible  to  drop  a perpendicular  from  any  point  on  the 
wire  to  the  plate  and  are  at  least  a quarter  wavelength  away  from  the  plate 
edges,  then  the  currents  may  be  found  quite  accurately  by  replacing  the  finite 
polygonal  plate  with  an  infinite  plate  and  then  solving  the  resulting  problem 
using  moment  methods  [l].  Many  computer  codes  exist  which  are  capable  of 
solving  such  a problem.  One  such  code  is  the  Antenna  Modeling  Program  (AMP) 
[2j.  The  input  format  for  our  program  conforms  to  that  of  the  output  from 
AMP.  It  is  assumed  that  the  wire  segments  are  small  enough  so  that  the  cur- 
rent flows  axially  along  the  wires  and  is  constant  over  the  length  of  one 
wire  segment. 

The  flat,  perfectly  conducting  polygonal  plate  may  have  up  to  50  i 

straight  edges.  The  interior  angles  made  by  the  intersecting  edges  must  be 
less  than  180°,  i.e.,  the  plate  must  be  a convex  polygon.  The  edges  can  be 
of  any  length  and  the  wedge  angle  of  each  edge  can  be  specif ed  by  the  user. 
This  enables  the  user  to  model  simple  solid  bodies. 

The  user  may  choose  not  to  include  the  diffraction  from  one  or  more  of 
the  edges.  This  option  makes  it  possible  to  model  semi-inf inite  structures, 
since  the  diffraction  from  the  edges  at  "infinity"  may  thus  be  excluded. 

The  infinite  surface  in  the  z = 0 plane  can  be  used  in  three  ways. 

First,  it  may  be  considered  as  an  impedance  surface  and  its  complex  normalized 
surface  impedance  specified.  The  normalized  surface  impedance  zg  is  defined 
by  Ey  = -zsn0Hx,  Ey  = zsh0Hy,  where  Ex,  E , Hx>  and  Hy  are  evaluated  on  the 
x-y  plane  and  nQ  = 377  ohms.  Special  provisions  for  specifying  the  surface  as 
a perfectly  conducting  groundplane  have  been  made  to  avoid  possible  numerical 
difficulties.  Second,  the  infinite  half-space  z < 0 may  be  filled  with  mate- 
rial of  arbitrary  permittivity  and  conductivity.  This  option  might  be  used 
in  modeling  structures  over  the  ocean  or  over  lossy  ground.  Note  that  for  any 
of  the  above  cases  only  the  fields  in  the  z > 0 half-space  are  determined. 
Thirdly,  the  surface  (and  the  entire  z < 0 half-space)  may  be  free  space. 

One  of  the  options  available  to  the  user  is  to  include  only  the 
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Figure  2. 


A square  plate  in  the  y-z  plane  illuminated  by  a current  element 
centered  over  the  plate.  The  frequency  is  10.43  gHz. 
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scattered  fields  when  summing  the  field  contributions  from  the  plate  and  sur- 
face. In  this  mode  one  (or  more)  current  elements  can  be  used  to  illuminate 
the  plate/ground  surface  with  the  desired  incident  field.  The  program  then 
calculates  the  scattered  fields.  Note  that  in  this  mode  of  operation  the 
simultaneous  use  of  a program  such  as  AMP  to  find  the  segment  currents  is  not 
required . 

Another  option  available  to  the  user  is  to  specify  calculation  of  only 
the  diffracted  fields.  By  comparing  the  diffracted  field  pattern  with  the 
total  field  pattern  for  the  same  geometry,  one  can  quickly  determine  the  dom- 
inant radiation  mechanism  for  various  pattern  regions. 

The  user  can  control  tne  output  of  the  program  by  means  of  various  input 
constants.  If  fields  in  dB  are  desired,  one  can  obtain  separately  either  the 
theta  component,  the  phi  component,  or  the  total  power  of  the  radiated  field 
in  a given  direction.  This  flexibility  allows  the  user  to  obtain  an  output 
which  is  not  cluttered  with  unwanted  data  when  only  one  polarization  is  of 
interest.  However,  if  desired,  all  of  the  above  may  be  printed  out.  The 
normalization  of  the  maximum  dB  level  may  also  be  controlled  by  means  of 
input  constants.  If  all  three  quantities  are  printed  out,  then  all  will  be 
normalized  in  the  same  way  so  that  comparisons  between  the  three  will  be 
meaningful . 

The  complex  electric  field  strengths  in  volts/meters  for  the  two  polar- 
izations (9  and  <)>)  can  also  be  printed  out  for  any  point  in  space,  provided 
the  point  is  in  the  far  field  of  the  radiating  structure. 

Provisions  have  also  been  made  for  punched  card  output  and  for  the  con- 
trol of  possible  visual  displays  of  the  calculated  patterns  which  may  be 
implemented  at  RADC. 

The  patterns  may  be  computed  with  a fixed  9 and  variable  $,  or  vice 
versa.  Additionally,  the  incremental  angle  between  pattern  data  points  may 
be  controlled  by  the  user. 

A simple  example  of  the  type  of  geometry  which  can  be  handled  by  this 
program  is  shown  in  Figure  2.  A radiating  current  element  is  located  above  a 
square  plate.  In  Figure  3,  the  calculated  radiation  pattern  in  the  x-y  plane 
is  compared  with  measured  data  [3j.  The  agreement  is  seen  to  be  quite  gooc . 


II.  PROGRAM  LIMITATIONS 

The  GTD  is  basically  a high  frequency  method.  The  user  should  be  aware 
of  the  basic  limitations  of  the  present  computer  code,  which  includes  only 
singly  diffracted  rays.  This  means  that  for  accurate  results  the  plate  edges 
must  be  larger  than  one  or  two  wavelengths  and  the  current  segments  should  be 
at  least  X/2  away  from  the  nearest  edge.  The  plate  size  could  be  signifi- 
cantly reduced  in  an  improved  version  including  doubly-diffracted  edge  and 
vertex  diffracted  rays. 
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The  omission  of  vertex  diffracted  rays  from  the  computer  code  does  not 
significantly  affect  pattern  calculations  such  as  the  one  in  Figure  3 where 
the  plane  of  the  pattern  cut  is  not  near  any  of  the  plate  corners.  However, 
in  Figure  4 we  show  a pattern  calculation  for  the  same  geometry  as  Figures  2 
and  3,  but  with  the  pattern  cut  in  the  plane  of  the  plate.  The  discontinu- 
ities in  the  calculated  patterns  are  caused  by  the  lack  of  an  adequate  mathe- 
matical model  for  the  diffraction  from  the  corners  of  the  plate.  Useful  data 
can  still  be  obtained,  but  the  user  of  the  program  should  expect  pattern  dis- 
continuities when  a diffraction  point  is  near  a plate  corner. 


metal-roofed  building. 


As  is  the  case  with  any  numerical  solution,  numerical  difficulties  can 
arise  if  some  of  the  distance  or  length  parameters  of  the  problem  vary  from 
one  another  by  many  orders  of  magnitude.  To  avoid  these  problems  the  user  is 
cautioned  to  keep  the  plate  and  sources  within  1,000  wavelengths  of  the  ori- 
gin and  to  let  the  plate  edges  and  source  distances  differ  by  no  more  than 
4 orders  of  magnitude. 

When  using  the  program  to  calculate  diffracted  or  scattered  fields  far 
from  the  plate,  one  may  find  a discontinuity  in  the  calculated  field  in  the 
direction  of  the  specular  reflection.  This  is  due  to  the  vanishingly 
small  angular  width  of  this  specular  contribution  in  the  far  field,  and  the 
calculated  field  at  that  one  angle  is  incorrect.  The  contribution  is  retained 
in  the  solution  because  it  is  needed  in  order  to  obtain  correct  results  closer 
to  the  scatterer. 

If  the  wire  radiator  is  located  between  the  plate  and  the  ground  plane 
(impedance  surface  or  dielectric  halfspace) , the  code  will  produce  a calcu- 
lated pattern;  however,  several  types  of  rays  which  are  not  significant  when 
the  source  is  above  the  plate  may  be  significant  when  the  source  is  between 
the  plate  and  the  ground  plane.  Examples  of  such  rays  include  rays  reflected 
by  the  ground  plane  and  diffracted  by  the  plate,  rays  multiply  reflected  by 
the  plate  and  the  ground  plane,  and  such  rays  eventually  diffracted  by  the 
plate.  Figure  1 illustrates  all  of  the  types  of  rays  presently  included  in 
the  code.  If  the  problem  geometry  is  such  that  other  types  of  rays  make  sig- 
nificant contributions  to  the  pattern,  then  the  results  obtained  using  this 
code  will  be  more  or  less  in  error,  depending  on  the  relative  importance  of 
the  the  missing  rays,  and  should  be  used  with  caution.  It  should  be  empha- 
sized here  that  for  the  cases  where  the  source  is  above  the  plate  (or  where 
no  z = 0 plane  discontinuity  is  present) , the  code  does  include  all  of  the 
significant  rays  (subject  to  the  constraints  on  plate  size  and  corner  dif- 
fraction previously  mentioned) . 

Despite  the  limitations  listed  above,  the  present  code  is  capable  of 
calculating  radiation  patterns  for  some  fairly  realistic  geometries.  Two 
examples  of  this  are  shown  in  Figures  5 and  6. 


III.  PROGRAM  DESCRIPTION 

A description  of  the  operation  of  the  present  computer  code  will  be 
given  in  this  section.  A top  level  flow  chart  for  the  program  is  presented, 
and  brief  functional  descriptions  of  the  important  subroutines  will  be  pro- 
vided . 


As  previously  indicated,  all  of  the  types  of  rays  included  in  the 
present  code  are  shown  in  Figure  1.  The  computer  code  thus  has  two  basic 
functions.  First,  it  must  decide  whether  or  not  any  ray  (or  rays)  of  a given 
type  exist.  If  so,  it  must  then  evaluate  the  electric  fields  associated  with 
these  rays.  The  results  of  the  above  process  for  each  current  segment  are 
then  added  to  produce  the  desired  radiation  pattern.  It  is  evident  that  the 
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computations  can  actually  be  separated  into  two  distinct  parts:  the  geometri- 
cal problem  of  finding  the  contributing  rays  and  the  electromagnetic  problem 
of  evaluating  the  associated  fields.  We  will  see  that  the  organization  of 
the  code  is  such  that  in  many  respects  the  two  parts  are  treated  distinctly, 
with  different  subroutines  being  used  for  the  geometrical  and  electromagnetic 
calculations . 

A top  level  flow  chart  of  the  main  program  (subroutines  are  excluded) 
is  shown  in  Figure  7.  From  this  it  is  evident  that  most  of  the  actual  calcu- 
lations are  performed  in  the  subroutines.  The  various  subroutines  are  ex- 
tensively documented  with  comments  which  are  distributed  throughout  the  pro- 
gram listing.  However,  in  order  to  provide  an  understanding  of  how  the  sub- 
routines are  logically  interconnected,  we  will  give  brief  descriptions  of  the 
basic  subroutines  ENGO,  ETDIF  and  GROUND,  list  the  subroutines  which  they 
call,  and  give  brief  descriptions  of  their  purpose.  These  descriptions  will 
be  given  in  brief,  blocked  out  form  for  easy  reference. 

ENGO:  Calculates  the  direct  and  incident  radiation 

fields. 

ENGO  calls  subroutines  EXIST,  IMAGE,  SOURCE. 

EXIST:  Determines  whether  or  not  the  incident  ray 

or  reflected  ray  exists. 

IMAGE:  Determines  the  location  of  the  image  of  a 

point  through  a plane. 

SOURCE:  Can  compute  either  1)  the  incident  field  at  a 

diffraction  point  or  2)  the  far  field  of  a 
source  or  image.  When  called  by  ENGO,  function 
2)  is  utilized. 

Subroutine  ENGO  loops  over  all  of  the  current  segments.  For  each  current  seg- 
ment it  cills  IMAGE  to  locate  the  image  of  the  source  in  the  plate,  then  calls 
EXIST  to  see  if  the  ray  from  this  image  passes  through  the  plate.  If  so, 
SOURCE  is  called  to  evaluate  the  field  from  the  image  location.  To  evaluate 
the  direct  radiation  ENGO  calls  EXIST  again  to  see  if  the  ray  from  the  source 
passes  through  the  plate.  If  not,  SOURCE  is  called  to  evaluate  the  field 
from  the  source  location.  The  two  contributions  are  added,  and  the  next  cur- 
rent segment  in  the  DO  loop  is  considered. 

ETDIF:  Calculates  the  total  diffracted  field. 

ETDIF  calls  subroutine  EMDIF. 

EMDIF : Calculates  the  field  diffracted  (if  any)  by 

a given  edge  in  a given  direction  when  illumi- 
nated by  a source  at  a given  location. 

Subroutine  ETDIF  loops  over  all  of  the  current  segments  and  all  of  the  dif- 
fracting plate  edges.  Each  time  EMDIF  is  called  to  evaluate  the  diffracted 
field  (if  any).  The  results  from  EMDIF  are  summed  to  give  the  total  diffrac- 
ted field.  (By  total  diffracted  field,  we  mean  the  total  field  diffracted 
from  the  edges  of  the  plate  to  the  field  point  due  to  direct  illumination  of 
the  edges  by  the  current  segments.  Doubly  diffracted  fields  or  fields  due  to 
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Figure  7.  Top  level  flow  chart  for  the  GTD-AMP  computer  code 
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i Figure  7 (continued).  Top  level  flow  chart  for  the  GTD-AMP 

computer  code. 
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reflections  from  the  impedance  surface  are  not  included) . 

GROUND:  Calculates  the  fields  reflected  from  the 

impedance  surface  whether  due  to  direct 
illumination  by  a current  segment  or  il- 
lumination due  to  diffraction  by  a plate 
edge . 

GROUND  calls  subroutine  RFPN,  SOURCE  and  EMDIF. 

RFPN : Determines  whether  or  not  the  ray  reflected 

by  the  impedance  surface  in  a given  direction 
can  exist. 

Subroutine  GROUND  loops  over  all  of  the  current  segment  sources.  It  calls 
subroutine  RFPN  to  determine  whether  the  reflected  ray  from  a particular 
source  exists.  If  so,  subroutine  SOURCE  is  called  to  evaluate  the  field 
from  that  particular  source.  Subroutine  GROUND  then  loops  over  all  of  the 
edges  (this  loop  is  nested  inside  the  source  loop) . For  each  edge  EMDIF  is 
called  to  calculate  the  dif f racted-ref lected  field.  Subroutine  EMDIF  itself 
makes  use  of  the  logic  in  RFPN  to  determine  whether  or  not  a given  diffracted- 
ref lected  ray  can  exist. 

The  above  information  should  provide  the  reader  with  an  understanding 
of  the  basic  organization  of  the  program.  For  more  detailed  information  the 
reader  is  referred  to  the  program  listing  itself.  As  previously  mentioned, 
this  listing  is  liberally  interlaced  with  comments  which  explain  the  working 
of  the  program  in  detail. 


IV.  INPUT/OUTPUT  DESCRIPTION 

The  various  coordinate  systems  which  are  used  to  specify  the  plate 
geometry  and  field  point  location  are  shown  in  Figure  8.  The  usual  Cartesian 
(xyz)  coordinates  are  used  to  locate  each  of  the  vertices  of  the  flat  poly- 
gonal plate.  Each  vertex  is  given  a number,  and  the  vertices  are  numbered 
consecutively  (in  either  direction)  around  the  plate.  Each  of  the  plate  edges 
is  also  assigned  a number  (shown  in  parenthesis) , with  edge  number  K between 
vertices  K and  K + 1. 

As  mentioned  previously,  the  diffracted  fields  from  individual  edges 
may  be  excluded  from  the  calculations.  This  makes  it  possible  to  treat  semi- 
infinite structures.  For  example,  the  model  of  an  infinitely  long  strip 
shown  in  Figure  9 has  the  length  of  the  strip  much  greater  than  the  distance 
from  the  source  to  the  strip,  and  the  program  can  be  directed  to  neglect  the 
diffraction  from  edges  (2)  and  (4) . It  should  be  noted  that  the  plate  must 
always  have  at  least  1 diffracting  edge. 

The  usual  spherical  coordinates  (r,  0,  if)  are  also  indicated  in  Figure 
8,  and  these  together  with  the  i (/  coordinate  are  used  to  specify  the  field 
point.  The  program  allows  for  two  types  of  pattern  cuts:  1)  conical  (or 

planar)  cuts  where  0 is  held  constant  with  varying  between  0°  and  360°,  and 
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— — ra&agr.: 


Flat  Plate 


/ 


180° 

360° 

360° 


Figure  8.  Coordinates  used  to  specify  the  plate  vertices, 
current  elements  and  field  point. 

2)  planar  cuts  with  i|>  varying  between  0°  and  360°.  The  additional  coordinate 
1 1>  is  introduced  to  avoid  the  difficulty  involved  in  making  a complete  circle 
around  the  radiating  structure  with  4>  constant.  This  difficulty  arises  be- 
cause 0 is  usually  constrained  so  that  0 - 0 - 180°.  The  plane  of  the  pattern 
cut  for  patterns  of  type  2)  is  determined  by  the  z axis  and  the  point  speci- 
fied by  the  r value  and  the  initial  and  4>  values.  This  point  also  deter- 
mines the  positive  direction  in  the  plane  (i.e.,  from  the  negative  z axis 
toward  the  initial  point).  The  value  of  r given  for  the  field  point  should 
be  in  the  far  zone,  i.e.,  it  should  be  far  enough  away  from  the  radiating 
structure  so  that  all  rays  from  the  structure  toward  the  field  point  are  n 
essentially  parallel.  Generally,  this  condition  can  be  obtained  if  r > 2D  /X, 
where  D is  the  maximum  extent  of  the  radiating  system  including  all  wire 
elements,  plate  edges,  and  all  image  points.  Phase  is  referenced  to  the 
origin  of  the  xyz  coordinates. 
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The  center  of  each  radiating  current  segment  is  specified  using 
Cartesian  coordinates  as  illustrated  in  Figure  10.  Either  the  unprimed  (xyz) 
coordinates  used  to  locate  the  plate  vertices  or  the  primed  (x'y'z*)  coordi- 
nates defined  by  the  plate  location  may  be  used.  This  second  coordinate 


X 


Figure  10.  Coordinates  used  to  specify  wire  segment  locations. 

Either  the  unprimed  (xyz)  coordinates  used  to  locate 
the  plate  or  the  primed  (x'y'z1)  coordinates  defined 
so  that  the  plate  is  in  the  x'y’  plane  may  be  used. 
The  x"y"z"  axes  used  to  define  angles  and  are 
parallel  to  whichever  coordinate  system  (primed  or 
unprimed)  is  chosen. 


system  is  added  so  that  coordinates  from  the  AMP  program,  which  requires 
that  the  ground  plane  coincide  with  the  x-y  plane,  can  be  fed  directly  into 
the  present  code.  The  primed  coordinates  are  defined  with  z*  perpendicular 
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to  the  plane  of  the  plate  and  x'  along  edge  (1),  with  the  origin  at  vertex  1. 

The  orientation  of  each  segment  is  specified  by  the  a and  6 angles 
shown  in  Figure  10,  where  x",  y"  and  z"  are  parallel  to  the  coordinate  system 
axes  used  to  specify  the  segment  center  (either  x,  y,  z or  x*,  y’,  z').  Hiis 
specification  of  the  wire  segment  location  is  the  same  as  that  used  in  the 
output  of  AMP. 

The  connection  data  and  tag  numbers  used  by  AMP  may  be  input  to  the 
code,  if  desired.  This  information  will  then  be  printed  out  for  identifica- 
tion purposes.  However,  this  data  is  not  needed  or  used  by  the  OSU  code. 

We  will  now  describe  in  detail  the  format  of  the  input  cards  which  con- 
trol the  operation  of  the  code.  The  input  data  will  consist  of  a section  of 
comments  and  four  numerical  input  sections.  The  comments  section  is  a set  of 
comment  cards  which  enable  the  user  to  print  out  a description  of  the  problem. 
Each  comment  card  starts  at  column  1 with  an  asterisk  sign  (*) , except  the 
last  comment  card  which  must  begin  with  a period  sign  (.).  At  least  one  com- 
ment card  must  be  included  in  the  input  data  set.  After  the  comments  are  the 
following  sections: 

-1  Plate  data, 

-2  Source  data, 

-3  Pattern  control, 

-4  Output  control. 

All  four  sections  need  not  be  included.  However,  those  included  must  be  in 
the  above  order. 

The  first  card  of  each  section  identifies  the  section  number  at  column 
15.  Columns  1-10  are  for  identification,  which  can  contain  any  combination 
of  characters  to  help  the  user  to  identify  the  card.  The  contents  of  columns 
1-10  do  not  have  any  effects  on  the  computer  code  and  will  be  referred  to  as 
the  identification  field.  A number  "1"  in  column  20,  whenever  it  appears  on 
the  first  card  in  a section,  indicates  that  the  program  will  be  executed 
right  after  reading  that  section's  data;  otherwise,  the  data  of  the  next  sec- 
tion will  be  read  following  that  of  the  present  section.  For  more  details, 
refer  to  the  discussion  on  modification  of  input  data. 

A.  Plate  Data 

Following  the  plate  data  section  card  (1  in  column  15)  is  the  plate 
data,  which  includes  three  types  of  cards  besides  the  section  card.  These 
three  types  of  cards  appear  in  the  following  order: 

— One  card  to  specify  the  number  of  edges  of  the  plate  (M) 
and  the  number  of  diffracting  edges  (ME) . It  is  usual  to 
have  M ■ ME,  except  for  some  special  cases  in  which  the 
diffracted  fields  from  certain  edges  are  not  included. 
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The  format  is  (lOx,  215);  lOx  is  the  field  used  for  identi- 
fying comments. 


— A second  card  (included  only  when  ME  is  less  than  M)  is  to 
specify  the  diffracting  edges  by  edge  numbers,  one  number 
for  each  diffracting  edge.  The  format  is  (lOx,  1415)  and 
it  can  be  continued  on  other  cards. 

— A third  type  of  card  (one  card  for  each  vertex)  is  to  specify 
the  positions  (x,  y,  z)  of  the  vertices  and  the  corresponding 
wedge  angle.  The  format  is  (lOx,  4F10.5).  The  position  of 
one  vertex  and  the  wedge  angle  for  the  edge  between  that  ver- 
tex and  the  next  is  contained  on  one  card.  The  ordering  of 
vertices  is  described  in  Figure  2 and  must  be  applied.  All 
dimensions  are  in  meters. 

— The  last  card  in  this  section  is  used  to  specify  the  impedance 
surface  parameters.  The  Format  is  (lOx,  15,  2E12.4).  The 
integer  variable  IGND  in  the  15  field  may  be  coded  in  one  of 
three  ways;  IGND  = 0,  no  impedance  surface,  i.e.,  free 
space  everywhere;  IGND  = 1,  dielectric  half space  for  z < 0; 
IGND  = 2,  impedance  surface  at  z = 0.  For  IGND  = 1,  the 
relative  permittivity  Er  and  the  conductivity  o (mhos/meter) , 
respectively,  are  coded  into  the  two  E12.4  fields.  For 
IGND  = 3,  the  real  and  imaginary  parts  of  the  surface  imped- 
ance zs  (normalized  to  that  of  free  space)  are  coded  into 
the  two  E12.4  fields.  However,  if  e and  o (or  z ) are 
equal  to  zero,  the  z = 0 surface  will  be  treated  as  a 
perfectly  conducting  ground  plane.  (These  zero  values  are 
then  not  related  to  physical  properties,  but  serve  only  to 
input  into  the  code  the  information  that  the  x-y  plane 
should  be  treated  as  a perfect  conductor) • 


B.  Source  Data 

Following  the  Source  Data  Section  card  (2  in  column  15)  are  three  types 
of  data  cards: 

— Number  of  current  segments  |n|  and  coordinate  system  used  — 

N negative  for  primed  system;  positive  for  unprimed.  Format 
(lOx,  15). 

— Segment  positions:  one  card  for  each  segment.  On  each  card: 

position  of  midpoint  (x,  y,  z)  or  (x1,  y',  z'),  length,  angles 
alpha  and  beta,  and  the  connection  data  (optional) . Format 
(6F10.5),  415). 

— Complex  segment  currents,  rectangular  form,  format  (2E12.4). 

Again  there  is  one  card  for  each  current  segment.  The  order 
must  correspond  with  that  of  the  above  segment  position  cards. 
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C.  Pattern  Control 


Following  the  Pattern  Control  Section  card  (3  in  column  15)  is  the  in- 
put card  which  controls  the  pattern  calculations  and  frequency.  It  is  for- 
matted (lOx,  5F10.5,  Ell. 4).  The  input  variables  and  descriptions,  in  the 
order  in  which  they  are  read,  are  listed  below: 


AI 

AF 

DEL 

PHITHE 


FR 


RA 


the  initial  ($  or  < (i)  pattern  angle 
the  final  (4>  or  i)>)  pattern  angle 
the  angle  increment 

the  fixed  (9  or  $)  pattern  angles.  If  PHITHE 
is  negative,  it  denotes  a fixed  0;  if  positive 
or  zero,  a fixed  <(>.  (0  = |PHITHE|). 

frequency  in  gigahertz;  if  FR  = 0,  the  default 
frequency  of  .3  gHZ  is  used  (the  wavelength  is 
then  1 m,  so  that  all  lengths  can  be  input  as 
wavelengths) . 

radial  distance  from  center  of  coordinate  system 
to  field  point;  if  RA  = 0,  then  the  default  value 
of  10,000  wavelengths  is  used.  Note  that  the  radial 
distance  specified  must  be  in  the  far  field  of  the 
source-plate  geometry  for  accurate  results. 


D.  Output  Control 

The  output  control  portion  of  the  input  data  consists  of  the  section 
card  (4  in  column  15)  followed  by  one  card,  formatted  (lOx,  815,  F10.5).  The 
control  variables  on  this  card,  in  the  order  in  which  they  are  coded,  are: 

— IOUT1 : If  only  the  0 component  pattern  or  only  the  4> 

component  pattern  or  only  the  power  pattern  is  desired, 
then  this  variable  is  set  to  a non-zero  integer  value  as 
follows: 

1 for  0 component  of  pattern  in  dB, 

2 for  $ component  of  pattern  in  dB, 

3 for  power  pattern  in  dB. 

The  individual  component  selected  will  have  a max  value  i 
in  dB  equal  to  FACTOR  (described  below) . If  IOUT1  = 0, 
none  of  the  individual  patterns  listed  above  are  printed 
out. 


— IOUT2:  If  all  three  of  the  above  patterns  are  of  interest, 

then  set  I0UT2  = 1.  This  will  result  in  the  0 and  $ com- 
ponent patterns  and  the  power  pattern  being  printed  out  in 
dB.  All  three  will  be  normalized  identically,  so  that  com- 
parisons between  the  three  components  will  be  meaninqful. 
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The  maximum  value  of  the  power  pattern  will  be  equal  to 
FACTOR.  If  IOUT2  = 0,  the  three  components  will  not  be 
printed  out. 


— I0UT3:  If  complex  electric  fields  are  desired,  then  set 

I0UT3  as  follows: 

I0UT3  = 1 will  output  the  complex  0 and  $ 
components  of  the  electric  field 
in  rectangular  form 

I0UT3  = 2 will  output  the  complex  0 and  $ 
components  of  the  electric  field 
in  polar  form 

If  the  complex  electric  fields  are  not  desired, 

set  I0UT3  = 0. 

— J0UT1 , J0UT2,  J0UT3 : If  any  of  these  variables  are  set  = 1, 

the  corresponding  (I0UT1,  I0UT2,  I0UT3)  field  quantities  will 
be  punched  on  cards.  If  card  output  is  not  desired,  set  them 
= 0. 

— J0UT4 : This  variable  reserved  for  control  of  plotting  routines 

which  may  be  implemented  at  RADC. 

— IFIELD : Determines  the  nature  of  the  calculated  fields: 

-1  for  diffracted  field  only, 

0 for  total  field, 

1 for  scattered  field  only. 

— FACTOR:  Maximum  value  for  dB  patterns  controlled  by  I0UT1 

and  I0UT2.  See  the  discussion  above. 


E.  Modification  of  Input  Data 

The  computer  code  is  designed  to  handle  many  different  problems  at  each 
run.  If  certain  data  sections  do  not  vary  from  one  problem  to  the  next,  the 
data  need  not  be  read  in  again.  After  the  first  problem,  the  data  set  need 
only  consist  of  the  introduction  (comments)  section  and  the  modified  data  sec- 
tions. The  last  section  of  modified  data  to  be  read  in  for  each  problem 
should  contain  a 1 punched  in  column  20  of  the  section  card.  The  remaining 
cards  of  that  section  will  be  read  and  the  program  will  then  begin  calcula- 
tions for  that  problem. 


V.  SAMPLE  COMPUTER  RUNS 

In  this  section  we  will  present  four  sample  runs  of  the  present  OSU 
program.  Each  sample  will  be  complete  with  input  cards  and  program  output 
listed.  The  first  three  samples  will  be  cases  where  the  desired  results  can 
be  obtained  without  using  AMP  to  find  the  currents.  The  last  case  will 
first  utilize  AMP  to  find  the  currents  (the  AMP  output  is  included)  and  then 
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the  currents  will  be  input  into  the  OSU  code  for  pattern  calculations. 


A.  Current  Element  Above  a Square  Plate 

The  first  sample  run  will  correspond  to  the  geometry  of  Figure  2,  and 
the  pattern  cut  will  produce  the  pattern  of  Figure  3.  The  input  card  listing 
and  program  output  are  shown  in  the  next  pages.  Since  this  is  our  first 
example,  a brief  discussion  of  the  input  data  stream  is  appropriate.  The 
prospective  user  should  follow  this  discussion  with  reference  to  the  previous 
section  of  this  report. 

The  first  five  cards  are  comment  cards  used  to  identify  the  run.  The 
last  comment  card  starts  with  a period.  Following  these  cards  the  actual  in- 
put is  divided  into  four  sections  by  the  section  cards  which  begin  with  the 
word  "SECTION"  (optional  — for  user  identification)  and  are  identified  by 
number  in  column  15. 

Section  1 is  concerned  with  the  plate  and  ground  plane  data.  For  this 
problem  the  plate  has  4 edges,  all  4 of  which  may  contribute  to  the  diffracted 
fields.  Each  vertex  is  given,  and  the  interior  wedge  angle  (left  blank)  is 
0.0°  for  each  edge.  No  impedance  surface  is  present. 

In  Section  2 the  number  of  current  segments  (1  in  this  case) , their 
position,  and  the  complex  current  flowing  on  each  is  given. 

In  Section  3 the  pattern  control  data  is  given.  For  this  run  0 is  fixed 
at  90°  (x-y  plane)  while  4>  varies  from  0°  to  360  in  5 increments.  The  fre- 
guency  and  the  distance  from  the  origin  to  the  field  point  are  given. 

In  Section  4 the  program  is  directed  to  produce  a printed  output  con- 
sisting of  the  0 and  $ components  in  dB  with  a maximum  value  of  00.0  dB  and 
also  the  complex  electric  field  (polar  form)  for  both  the  0 and  components. 
Hie  total  field  will  be  calculated. 

Ttie  corresponding  output  is  listed  following  the  input  card  listing. 


Case  A Input 


•PROBLEM  : 

•SHOOT  DIPOLE  ABOVE  FIAT  PLATE 
•PLATE  IS  I?"  SQUARE 

•DIPOLE  IS  6.75"  ABOVE  THE  PLATE  C IN  THE  CENTER 

.FREQUENCY  IS  10. A3  GHZ 

SECTION  1 

NO.  FOGES  A A 


VERTEX 

0.0 

0.15240 

-0.1524 

VERTEX 

0.0 

0. 15240 

0. 1524 

VERTEX 

0.0 

-.15240 

0.1524 

VERTEX 
NO  GROUND 

0.0 

0 

-.15240 

-0.1524 

SECTION 

2 

• SEGMENT 
0. 1714 

0.0 

1 

0.0 

0.01 

90. 

0. 

l.OOOOE+OO 

0. 

OOOOE  *00 

SECTION 

PATTERN 

O 

• 

o 

3 

360. 

5. 

1 

X> 

o 

• 

10.43 

SECTION 

4 

OUTPUT 

0 

1 2 

0 0 

0 0 

00.0 

♦1.0000E+05 
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Output 


I 


I 


!♦  AND  I-  INDICATE  THE  SEGMENTS  BEFORE  ANO  AFTER  I 
COORDINATES  OF  SEG.  CENTER  SEG.  ORIENTATION  ANGLES  CONNECTICN  DATA  TAG  * - CURRENT 


tmpta  i phi  i NO«*Ai!zPn  ppw«r  gain 

1 I E-THFTA  | F-PHI  | TPTftt 


NMT 
*■  r a 
f*  I Mf 


r\r*'  r~  o cr  — r-  t~  r-  4 r v _.f  r er  a ' r~  c o f o c +•  r o r -r  •♦  tr  r — r»i 

■C  T O'  »••  —«<r  oc  or  • Mfff(jM7irj>j®MOO  — J — •*-  — .*  r — r tr  or  r-  «■#  ■#> 

rrr  r-  — r r«r#r.j<r#»j.n  \ jr-  — r J «\<c  jcaa^rr-CaC'#-  ry.*-  — a rs.  ^C' 

•c.  a * r^o  a.  O ^ * r r t'  r*  ir  a r t <**  9 rv  rv-ir  at  a r-  J O J O Jr-  'r  a »o  r,  »»  r_  <y  irf  fv 


m • 1 i T i W i i T • • *i  *1 T TfT T i i i T i _r\*>  1 T i 1 1 f 

i i * 1 1 1 • i 1 1 1 1 1 * i * i 


r^fNiCff  -'--rrrf'r  X -c  ^ cor-  ► c*—  -c  -c  »r  -r  or  r-  — coNitm 

r- —• to  c,f\*roc- -cc  o Cc  <jo  * • — »_  »r  o <*  c j or-  r>o  a.  t j «t  a <■  «#  •# oo  4.  ir 

49'ru4«^M'fi<tirsf*>ai<4^u>^4i/'irirer/)ff  «r  xur  tor-'O-r  rJ»/>irrror'»»  jbc 

■for  ■£  —•€  a *'r\-% r tr  t o tr*  c o *r.r-t\C  c-  tr-  o c or  *rr*.<r<r  -rt^^r-  o o-~o  t r * c -tr- C r r-< 


—•to  f O -To  O -trj^r~~*r*'*cor<.-Crr  ot^  or s.  «r  *r  »*>  *-  P-  Jf\»  J c*  CffOf  *0  cr-  rv«  -*t*\  cr  or  «r«“  ffiMTOO1 
J j .0 .0  ■*.  r-  r~  o.r-  a r-r-  OOO  Of  Ou  tf'fs.i*  **>  OoO  00  X 

»-«.  « — « ■ iw^«^<Nf^»inf^«>.rgr^i  ^rvNi^nrvNNNrviNr^Nrwi  «,-»-  • -■»»-■  ^rwivi>iNN'Wi  i> 

t I I I I I I I ••  I I • I I I I I I t • I I I I I I I • I I I • I I I I I I I I • I I t I I I I I I 


r-  -cr*  jsf*o  -f  r ^rir  — <r  a 4^  .-*c  f tr  r~  c u*  <xrc  ^ ror  •— r ,4*r  o— «r- 

r^rvcr  ^orro-r^^ac  >0  acr-*<'0<?r-o‘r'jj  ccMroa-*-!  —r-c.  10-- a: -•iooic-~  j •— ao  im44< 
*-  f'a  *»M\ir--*c  r iTr  4arr'4<c  ^41441.  ~-«rv  <r  -.f-ror  -t  occ  or-  rrr-  o * 0*  on»>-« « rv— o ■© 
t\jt.r*c  «rr  M>a  4)^r  i»Ntnro’r  r r*“  c*  rvf^rs  ithtkk  jo  fp  oc  nrr.^  r p ^ n« 


ore-  O rr*-.  eg  r-  — «•/'  f*' j)  J tT'J  W'»r>x:«/'P'r^  .0  O rg  to  r\.  C"  f^r-  f- r-  onkm^(  cotrcjO  O fr cmO 

III  I | I I I I I I I I I I I I | I I I I I I I • I | I | I I t 

I I I I I I I I I I I I I I I I I 


c*«  iHVpOocc'Crr  c or c c 00  00 c rrcc  ronr  cc  CjC  r.oc  c>c  oor  c rrr  r cxj o 
c (.ouococc «.  c c»r,r~»ot  cu  c.cc.c  c ccc  rcccc  r ccr  ccrotot  cc  cc  oeoo* 

ntrr  rrir  rtrr/,',ir'~'irr  irrrf'f  irr  10c  irrtrt  trrirr^irr^ui^  i/'rv'''if  oko  c or  a e irrirri 
O O'  OC  —•—rvj«\,r'  rr  j -•  «r  »r  oo*-*~*«cO  O t C —*—.cvC\i»*>r“i  j J ir. 

— — — — "N/<\  ry<r^rsir^ (vr^r^rkg  rs#< 


cc.cc  oc  coa  ccc  c cccoc  occc  c.  ;cor  c*c  c c c cc  ici  C/crCccuti  c c 0* 
rrrcc  c 1 < '“'j1  » * r r r c occcccct  rr  » cr  rcrrcrc^cCcctcf  c.  c <*  r- 

r.c'r  o<  ci'rrtrt » r'c  <*rO(»  rcrcrffrccc  trr  r crr<  cm  « rocor  '*r  cro< 

ppcpap0ffpop^poppaappppopa»ffp»aopepffpepacpoe»a(raoc»»' 

o 

o 
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W 


m 


i 


r tcrr  t i rr-t  r^r-  r X--  O f*  a •xu'  i*' 

f n C r ■•  i cr~(  it  — f*  7 <*“ « X C ferv. 

,*  r r (••«  -T  'r'  u"  < »■  -►•rai*-  Of  <- 

xi/urr-i  — a rvi  rw 


o«r  -4  ir  X -z  »*'«/■  *r  o x X'»-,*‘-  f\»«-*i*  oo»^n- cc 

I I I I I l l I I I t I I I I • I I ill 


nc  p-  «\|f'  orrsi^i/'^  <r  ■ p-x  <*•*.«><  'for 

Xvf  P-  rf  1/  « a a *T  CC  0 f-  f\-  X X -4  M X Xf'ilf 
a-  7 (7  O rr«r  k a X r*  m x P*  Y nr  • X 


JCaOPCCf  ffffC  * Ya  mix  — 

XX  .4,  XX  ■«.•*.  P-  X X X-CxXiCOoi^xXlf 

rs,r^<<4  '\n-*r*  'v.NNt  -.Mr^r^rN.  .-4  ^ 

• I I I I I I I I I I I I I I t I I I t I I 
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B.  Current  Element  On  Step 


The  geometry  for  this  problem  is  shown  in  the  insert  of  Figure  11.  We 
wish  to  find  the  radiated  fields.  This  is  a two-dimensional  problem  which 
can  be  solved  quite  easily  using  the  GTD.  The  step  structure  can  be  modeled 
by  the  GTD-AMP  computer  code  by  letting  the  plate  form  the  top  surface  of  the 
step  and  the  infinite  ground  plane  form  the  lower  surface.  The  vertical 
surface  can  not  actually  be  modeled.  However^  all  of  its  effects  can  be  in- 
cluded by  specifying  that  the  plate  have  a 90  wedge  angle.  Since  the  origi- 
nal problem  has  only  one  diffracting  edge,  the  GTD-AMP  code  is  directed  to 
include  diffraction  effects  from  only  the  corresponding  plate  edge.  The 
other  edges  are  specified  to  be  far  from  the  current  element  so  that  the  semi- 
infinite nature  of  the  original  problem  is  closely  approximated.  The  input 
data  and  corresponding  output  appear  on  the  following  pages. 
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C.  3/4 X Vertical  Monopole  Over  a Dielectric 


The  geometry  of  this  problem  is  indicated  in  Figure  12.  The  current  on 
the  monopole  is  assumed  to  be  sinusoidal  so  that  AMP  is  not  needed  to  deter- 
mine the  currents.  The  monopole  is  divided  into  8 current  segments.  The 
monopole  is  over  a dielectric  half  space  with  er  = 4,  and  no  diffracting 
plate  is  present  in  the  original  problem.  However,  to  fulfill  the  program 
logic  requirements  that  a plate  with  at  least  one  diffracting  edge  must  be 
included  in  each  problem  a very  small  plate  far  from  the  monopole  is  included, 
and  this  plate  is  further  constrained  to  have  only  one  diffracting  edge  so 
that  its  effect  on  the  pattern  is  negligible.  The  results  from  the  GTD-AMP 
code  are  compared  with  results  from  Collin  and  Zucker  [4]  in  Figure  12. 
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Case  C Input 
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SECTION  9 

OUTPUT  0 1 2 0 0 0 0 00.0 
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Case  C Output 
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D.  Three  Element  Yagi  Over  a Finite  Ground  Plane  -- 
AMP  Interface 


In  this  section  we  are  concerned  with  finding  the  pattern  of  a three- 
element  Yagi  antenna  in  the  presence  of  a 13.X  square  conducting  plate.  The 
dimensions  of  the  Yagi  are  shown  in  Figure  13. 

In  order  to  find  the  currents  on  the  Yagi  elements,  the  AMP  program  is 

used  and  the  currents  are  found  with  the  Yagi  located  above  an  infinite 
ground  plane.  As  discussed  previously,  these  currents  will  be  affected  very 
little  by  the  presence  of  the  plate  edges  provided  the  edges  are  at  least 
A/4  away,  as  is  the  case  here.  The  AMP  program  output  [5]  is  listed  on  the 
following  pages. 

The  currents  from  the  AMP  program  output  were  input  to  the  GTD-AMP  code 
via  the  input  cards  listed  following  the  AMP  output  listing.  Two  different 
geometries  were  included:  1)  Yagi  over  13A  square  plate;  2)  Yagi  over  infin- 
ite groundplane.  Note  that  the  Yagi  currents,  since  they  are  the  same  for 

both  geometries,  need  not  be  repeated.  The  results  from  the  GTD-AMP  code  are 
listed  and  are  compared  with  one  another  and  with  the  pattern  calculated  by 
AMP  in  Figure  13.  Measured  data  for  this  problem  is  given  by  Kraus  [6],  and 
shows  good  agreement  with  the  calculated  values  in  Figure  13. 
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a OSU  - 13  X Square  Plate  X R 
R = I05X 


0.58  X 050X 


0.45  X 


Wire  Radius 
0.001  X 


Figure  13.  Calculated  patterns  for  a 3 element  beam  antenna 
over  a ground  plane  and  over  a finite  plate. 
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Case  D AMP  Output 


:»n  >»d  »cm»» 


Case  D Input 


•CURRENT  SOURCES  OVER  A PLATE 
.CURRENT  DISTRIBUTION  OBTAINED  FROM  AMP 
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APPENDIX  A 


MODELING  THE  SPECULAR  REFLECTION 
FROM  THE  GROUND  OR  OCEAN 


The  reflection  of  EM  waves  from  the  rough  ground  or  ocean  is  determined 
primarily  by  two  quantities  — the  constitutive  parameters  of  the  medium  and 
the  surface  roughness.  If  the  surface  is  relatively  smooth,  then  the  reflec- 
tion will  be  primarily  specular  and  can  be  determined  quite  easily  from  the 
complex  permittivity  of  the  medium.  If  the  surface  is  very  rough,  then  the 
reflection  will  be  incoherent  and  the  simplest  approximation  to  the  specular 
reflection  is  to  neglect  it  altogether.  If  the  surface  is  neither  very 
rough  nor  very  smooth,  the  problem  of  determining  the  reflection  is  quite  in- 
volved and  is  beyond  the  scope  of  the  present  work.  However,  for  certain 
cases  (for  example,  the  ocean)  an  effective  surface  impedance  can  be  defined 
which  will  allow  the  calculation  of  an  approximate  average  reflection  from 
the  rough  surface. 

A measure  of  surface  roughness,  which  can  be  used  to  quantitatively 
determine  relative  roughness  or  smoothness,  is  the  Rayleigh  criteria: 


4tis  cos0 


where  s = Standard  deviation  of  surface  irregularities 

0 = Angle  of  incidence  measured  from  the  normal 
\ = Wavelength. 

For  R < 0.1  there  is  a well-defined  specular  reflection  and  the  reflection  co- 
efficient may  be  approximated  using  the  usual  plane  wave  reflection  coeffi- 
cients (in  our  code,  this  corresponds  to  option  IGND  = 1 in  the  Plate  Data 
Section).  For  0.1  < R < 10  the  surface  roughness  may  not  be  neglected.  How- 
ever, if  an  effective  surface  impedance  is  available,  it  may  be  used  to  deter- 
mine an  approximate  reflection  coefficient  (option  IGND  = 2) . When  R > 10, 
the  surface  is  so  rough  that  the  specularly  reflected  field  is  negligible  and 
can  be  neglected  (option  IGND  = 0) . 

For  making  calculations  for  a "typical"  lossy  earth,  the  relative  di- 
electric constant  will  vary  from  about  7 for  a low  conductivity  earth  to 
about  30  for  a highly  conductive  earth,  with  a typical  value  being  er  = 15. 

The  corresponding  range  for  the  conductivity  is  from  1 to  30  mi 11 imho/meter 
[7]. 

For  frequencies  up  to  .5  gHZ  Barrick  (8]  has  developed  an  effective  sur- 
face impedance  which  can  be  used  to  approximate  the  reflection  from  the  ocean's 
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surface.  A graph  giving  values  of  this  normalized  surface  impedance  for  vari- 
ous states  of  sea  roughness  is  shown  in  Figure  A-l. 

Finally  it  should  be  emphasized  that  the  determination  of  the  reflection 
from  a statistically  rough  surface  (R  > 0.1)  is  a complex  problem,  so  that  the 
methods  outlined  here  for  dealing  with  this  problem  are  quite  simplistic  and 
can  be  expected  to  yield  only  approximate  results  at  best. 
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§ Rome  Air  Development  Center  ? 


RA DC  is  the  principal  AFSC  organization  charqed  with 
planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  information  sciences,  intelli- 
gence, command,  control  and  communications  technology , 
products  and  services  oriented  to  the  needs  of  the  USAF. 
Primary  RADC  mission  areas  arc  communications , electro- 
magnetic guidance  and  control,  surveillance  of  ground 
and  aerospace  objects , intelligence  data  collection  and 
handling,  information  system  technology , and  electronic 
reliability,  maintainability  and  compatibility . RADC 
has  mission  responsibility  as  assigned  by  AFSC  for  de- 
monstration and  acquisition  of  selected  subsystems  and 
systems  in  the  intelligence,  mapping,  charting , command, 
control  and  cormuni cat  ions  areas. 
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